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Abstract: We study the formation of marginally trapped surfaces in the head-on 
collision of two shock waves in de Sitter space-time as a function of the cosmological 
constant and the shock wave energy. We search for a marginally trapped surface on 
the past light cone of the collision plane. For space-time dimensions D > 3, there 
exists a critical value of the shock wave energy above which there is no trapped 
surface of this type. For D > 3, the critical value of the shock wave energy depends 
on the de Sitter radius, and there is no this type trapped surface formation for a 
large cosmological constant. For D = 3, the critical value of the shock wave energy 
is independent of the cosmological constant. At the critical point, the trapped surface 
is finite. Below the critical energy value, the area of the trapped surface depends on 
the cosmological constant and the shock wave energy. 
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1. Introduction 



It is a well-established experimental fact that our current universe is expanding with 
constant acceleration, which is well described by the extremely small cosmological 
constant A = 10 _47 GeV 4 , = 0.726 ± 0.015 PJ. There is a common opinion 
that two ultrarelativistic point particles in the Minkowski space-time can produce 
black holes |2|, [3|, This longstanding question has a purely theoretical meaning 
as well as an astrophysical one ||. In the framework of TeV-gravity ||, black hole 
production in collisions of particles with the center-of-mass energy of a few TeV and 
their experimental signatures |7j became the subject of numerous investigations 0, 
H [U], |ll|, [12| . We also note a discussion of the possible production of wormholes 
and other more exotic objects at the LHC JT^, [14|, |15| (see for a consideration of 
wormholes in astrophysics). 

Our main aim in this paper is to determine if the presence of a positive cosmolog- 
ical constant can influence the processes of black hole formation in a collision of two 
ultrarelativistic point particles. Intuitively, it is almost clear that a small cosmolog- 
ical constant cannot have a detectible influence. Meanwhile, a positive cosmological 
constant generates a repulsion of matter, and a critical value of the cosmological con- 
stant A should exist above which black hole formation is suppressed. Finding this 
critical value is worthwhile. We note that under the assumption of an asymptotically 
flat space-time, the presence of a trapped surface usually guarantees the existence of 
the event horizon |T7|, [18], 15, ^0| . 



We study the formation of marginally trapped surfaces in a head-on collision of 
two shock waves in the de Sitter space-time as a function of the cosmological constant 
and the square of the shock wave energy. For D = 4, we show explicitly that there 
exists a critical value for the ratio of the shock wave energy p and the cosmological 
radius 

above which there is no solution of the trapped surface equation. Here G is the 
gravitational constant and the cosmological radius a is related to the cosmological 
constant by a 2 = 3/A. Similar results were found for D > 4. There is a critical value 
for the shock wave energy above which there are no solutions of the trapped surface 
equation. This value depends on the cosmological radius. At the critical point, the 
trapped surface is finite. 

This effect is similar to the emergence of critical behavior with respect to the 
wave width in the transverse space when a marginally trapped surface is formed in 
the head-on collision of two shock waves in the Minkowski space-time ^Fl, This 



effect depends on the number of dimensions and is a Choptuik-like critical effect, i.e., 
it is similar to the critical Choptuik behavior in gravitational collapse |23] (see 
and the references therein). 
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The formation of a marginally trapped surface in the collision of gravitational 
shock waves in AdSo was recently studied in |25|, ^6| . These studies are aimed at 
better understanding the entropy production in relativistic heavy ion collisions due 
to black hole production in a dual description. Despite the absence of a holographic 
dual description of QCD, describing the colliding heavy ions in terms of colliding 
gravitational shock waves in the anti-de Sitter space-time was suggested [[Z7], 
Black hole formation in collisions of the dual of the nuclei in the bulk is interpreted 
as formation of a quark-gluon plasma |2{| |3lf| . In AdS, a dimension-dependent 
critical behavior with respect to the wave width in the transverse space in the for- 
mation of a marginally trapped surface in the head-on collision of two shock waves 
was recently found [pi]] . For D = 4 and D = 5, there exists a critical value of this 
width above which the trapped surface never forms. We note that in the AdS space- 
time, the obtained results are qualitatively the same as those obtained in the flat 
space-time. 

This paper is organized as follows. We start with the setup and recall some basic 
facts about the generalization of the Aichelburg-Sexl shock wave geometry [ 



33, EE I 



to nonexpanding shock waves propagating in D-dimensional space-times with the 
cosmological constant |35], [36] [37], BR [39] GO. Ell. E3L fj2|. In Sect. 3, we calculate the 



critical value of cosmological constant depending on the shock wave energy below 
which the trapped surface occurs and give the area of the trapped surface. 



2. Setup 

Our main aim in this section is to present the setup for studying the formation of 
closed trapped surfaces in the head-on collision of two shock waves in dS space (this 
may be compared with the setup used to study the case without a cosmological 



constant in [TT[] and the case with the negative cosmological constant in [GJ| ^5] |2T|). 
We briefly recall the results in |35, B9|, G3] for the geometry of a shock wave 



propagating in the D-dimensional dS space-time (see Appendix A). In terms of the 



dependent plane coordinates, (u,v,x), x 



[x 



x 



D 



), satisfying — 2uv + x 2 



(a is the cosmological radius), the line element of the shock wave space-time is 

ds 2 = -Idudv + dx 2 + F{x)5{u)du 2 . {2.1) 



The shock wave shape function F is a fundamental solution of the equation 

D-2" 



AgD-2 



F = -16V2nG D p5( 



n,n ) 



[2.2) 



where A s d-2 is the Laplace-Beltrami operator on a (D— 2)-dimensional sphere E> D ~ 2 , 
n = x/\x\, no is the location of the particle on the sphere, p is the energy of the shock 
wave, and Go is the D-dimensional gravitational constant. The a/2 in the right- 
hand side results from our choice of the coordinate system (if we impose —dudv 
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instead —Idudv in the expression for the linear element, then it disappears). This 
metric is a solution of the Einstein equations for an energy-momentum tensor with 
the single nonvanishing component T uu ~ p5(u). In the standard parameterization 
of the (D— 2)-dimensional sphere by spherical angles . . . , i?d-2> the shock wave 
shape function corresponding to an ultrarelativistic point particle depends on only 
one spherical angle $d-2- The operator A s d-2 acts on F = F($ D _ 2 ) as 



di} D -2 dv D - 2 



A gD - 2 F = ^ sin 3 - D (^ D _ 2 ) ( — sin^-^^.^-j^ ) F(tf D _ 2 ) (2.3) 



see 



37j regarding shock waves in dS/AdS with multipole structures). 
For D = 4, we deal with 



(2.4) 



where £ = x 4 /a = cos $2 and p = pG^ is the rescaled energy. 

We now consider a collision of two waves of the type described above. We 
suppose that in the region {u < 0} U {v < 0}, i.e., the part of the space-time before 
the collision, the metric is given by 

ds 2 = -2du dv + dx 2 + F(£, ^)5(u)du 2 + £ 2 )8{v)dv 2 , (2.5) 

where £1 and £ 2 are the locations of the two colliding particles (see for the explicit 
formula for F(£,£j)). In independent coordinates (see Appendix C.I), the metric is 

2 -2dw da + 2d( d( + 2ff 1 (C, C) S(w)dw 2 + 2H 2 ((, Q 5{a)da 2 

[1 - (wa - CC)/2aT ' { ] 

where F;(C, () = H((, (, Q, Q, d = C(&> &), and 

*(C, f , 0, 0) = ff(cf) = i (1 + ^cc) f (l^fg) • (") 

A rigorous analysis of the formation of black holes in collisions would require 



solving the Einstein equations in the interaction region {w > 0, a > 0} (see, e.g., ||7, 
^| and the references therein). 

A sufficient condition for a black hole to form in the asymptotically flat case 
is the existence of a marginally closed trapped surface at the hypersurface {w < 
0, a = 0} U {w = 0, a < 0} @, gj| [TJ, g|, |||. We note that in non-asymptotically 



flat cases, there are no general theorems, but there is a common opinion that the 
existence of a marginally trapped surface can be used as an indication of black hole 
formation. 

In the coordinates used in line element (|2.6|) null geodesies are discontinuous 
across the wave fronts, w = and a = (see Appendix C) . So using such coor- 
dinates to find this trapped surface equation is inconvenient and can be avoided by 
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switching to a new coordinate system (W, E, T, T) in which the delta function terms 
are eliminated in the metric and the geodesies are continuous. Similar to the D'Eath 



and Payne coordinates P5 |, which are closely related to the explicit form of geodesies 
in the Minkowski space-time with the shock wave Ell. |33| , these coordinates are also 



closely related to the geodesies. This is a reason for us to study geodesies in the dS 
space-time with the shock wave (see Appendix B for details and references). 

In these coordinates the trapped surface that we seek has two parts, which are 
denoted here by Si and <S 2 in the respective regions £ < and W < 0. They are 
defined in terms of the two functions ^i(T, T) and ^(T, Y) by 



Si: 



w = 0, 

E = -* 1 (T,T), 



<S 2 



£ = 0, 



:2.8i 



with the additional boundary conditions at the shock wave intersection C C {W 
E = 0} 



*i(T,T) =0, * 2 (T,T) 



0. 



and 



1. 



(2.9) 



(2.10) 



Because Si and S2 are in the respective regions W < and E < 0, we also have 
*i(T,T) > and * 2 (T,T) > 0. 

The two functions ^ r i(T, T) and ^(Y, T) must be determined by imposing the 
condition that the surface they define is marginally trapped |18], |48|, [46fl , i.e., that the 
congruence of outgoing null geodesies orthogonal to the surface has zero expansion. 

The zero convergence equation for the D=4 case reduces to the equation (see 
Appendix D) 

+ 4W(T, T) = 0, (2.11) 



A, 



where 



A, 



2 1 + 



2^ 



dyd- 



(2.12) 



is the Laplace-Beltrami operator in complex coordinates (see Appendix D) and the 
functions 4>i,2 are related to the functions 2 defining the two halves of the trapped 
surface and the shock wave shape functions Hi 2 (cf. similar equations in the AdS 
?Tfl) 

2*i )2 - #1,2 



case 



25 



h > 2 = ITrW (2 - 13) 

The next question to be discussed is whether a black hole forms as a result of 
the head-on collision of two waves of the type described above. Head-on collisions 
preserve rotational symmetry around the axis of motion of massless particles, the 
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0(2)-symmetry in D = 4. Because a head-on collision is 0(2)-symmetric, the func- 
tions ^i(T, T) and \I/2 (T , T) describing the trapped surface are identical and depend 
on only the parameter p 2 = YY: \&i(T, T) = ^(Y, T) = ^(p 2 ). 



3. Marginally trapped surface in dS 



3.1 Solution of the marginally trapped surface equation in dS4 

In the case </>i(T,T) = </> 2 (Y,Y) = </>(p 2 ), equation (|2.11| ) is transformed into the 
ordinary differential equation 



2a 2 



+ + -7T = 0. 



(3.1) 



dp 2 p dp 

The solutions of this equation are 

_ A(p 2 - 2a 2 ) + B((p 2 - 2a 2 ) hip + 4a 2 ) 
<j)= p 2 + 2a 2 

where A and B are constants. Because we are interested in a solution of a homogenous 
equation without singularities, we take B = and obtain the expression for the 
trapped surface function in terms of the shape function H, 



(3.2) 



1 



V = lH(p) + ( 1 + 



2a' 



A{p 2 - 2a 2 ) 



and also in terms of the shape function F, 
or, more explicitly, 



2a 2 



(3.3) 



(3.4) 



(3.5) 



As mentioned above, the function \l/ must satisfy the following boundary condi- 
tions in the head-on case: 



d r ^df * 



0. 



1. 



(3.6) 
(3.7) 



It is obvious that the bound C is a circle p = po = const. We hence have a system of 
two equations for the two constants A and p : 



2a 2 



pi 



2a 2 + p 2 



2a 2 



2a 2 



0. 



^ ^v^pa 2 + p 2 (V2p -A) + v^ppg In ( 



2a 2 



(3.8) 
(3.9) 
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Substituting A from (3.8) in expression (3J3), we obtain 

-,2 



*(p) = ^2p (V* /; " 



V Po - 2a 2 



1 



P 
2a 2 



(3.10) 



Equation ( |3.9p defining p can be also written in terms of the initial function F 
defining the shock wave in the plane coordinates, 



A 

2a? 



F'(po) + 



Po 



2a 2 - p 2 



F(p ) + 2 = 0. 



(3.11) 



Indeed, normalization condition ( p,10| ) can be written in the form 

— = 4 or — = *'p =±2- 
\dp J dp 



(3.12) 



We choose the minus sign. Because of relation ( |3.4j ) on C, we have 



c 4 V" ' 2a 2 
and we obtain ( 2.10|) in the form 



Po 



2a 2 - p 2 



F(po 



(3.13) 



PL 
2a 2 



F'(po) 



Po 



2a 2 - p 2 



F(p ) + 2 = 0. 



(3.14) 



This equation is universal for an arbitrary dimension D. To connect with the AdS 
case, this equation can be rewritten in terms of the chordal coordinate related to p 
by 



p = a 



1-q 



After this change of variable, equation ( 3.1 1|) becomes 



F'(q ) + 



-F(q ) + 



0. 



(3.15) 



(3.16) 



l-2g " ^ ' ^2g (l - q ) 
We introduce the dimensionless parameter xq = Po/a. Equation (|3.9|) becomes 



f(xo) = V2- 



where 



P 



(3.17) 



(3.18) 



x 2 — x 2 

We note that in the region < x < \^2, the function f(x) has the positive 
minimum 

f(x) \ x=Xmin = 0, x min = 2-V2. (3.19) 
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Hence, in the case 

rj = - < -1= ■ f(x min ) = 4, (3.20) 
V V2 

there are no solutions of ( |3.i7| ) (see Fig. 1). Therefore, no solution of the trapped 
surface equation can be found below the critical value p c = 4 of the parameter p. In 
terms of the rescaling energy p = p/G, where G is the gravitational constant, this 
relation gives Ql.l ). 



It is reasonable to consider two limit cases: where xq <C 1 and p < a (when the 
particle energy is low and/or the spacetime is weakly curved; it is natural to call this 
the low-energy limit) and where Xq and rj are equal to the critical values (we call this 
the high-energy critical limit). 

From equation ( |3.17| ) in the region xq -C 1, we obtain 

\/2v 

Xq « and p ~ V2p. (3.21) 

a 

And in the critical limit, we obtain 

p = ~ Po = (2 - V2)a = (8 - ±V2)p. (3.22) 




Figure 1: The critical value of the parameter rj = a/p corresponds to the minimum value 
of the function / (red line). 



3.2 Higher-dimension cases 

Similar calculations can be performed in higher dimensions. Equation (|3.14|) has a 
universal form and is independent of the dimension D. The variable p for an arbitrary 
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dimension is related to £ = Zd /a as 

l-p 2 /2a 2 , , 

^TTTW- (3 ' 23) 

As examples, we consider the cases D = 5 and D = 6. The shape functions have 
the forms 

= 3^2^ (3 . 24 ) 

~p(wj + "^a" J ' (3 - 25) 



where £ = Z D /a and = pGp- For these cases, equation Q3.14Q reduces to 

foM = C D — , (3.26) 
Pd 

where Co is a D-dependent constant (in particular, C 3 = 2/tt, C 5 = 32/37T, and 
Cq = 6v2) and the functions /d(x) are given by 

These functions have positive minima at the points 



'2(-l + £> -2V-D -2) 

aw^D = — — , -D > 3, (3.28) 

V -L> — 3 

x miri , 3 = 0, -D = 3, (3.29) 

and the value of the functions /z?(x) at these points, / ,d = fD(%min,D), give the 
critical values of the parameter p below which trapped surfaces can be formed, 

P < Pcr,D, (3.30) 

a LjD lo ot\ 

Pcr,D = n , ■ (3.31) 

J0,D 

We note that this critical effect does not necessarily exclude the trapped surface in 
the region of interacting shock waves W > 0, S > 0. 

We can also interpret formula ( |3.31| ) to mean that for a fixed value of p, there is 
a critical value of the de Sitter radius a cr , 

r \ l/D-3 

Jo,d 



a cr = ( PGd-^ ) , D > 3, (3.32) 
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only above which can trapped surfaces be formed. For D = 3, the critical energy is 
independent of a. 

In particular, for D = 5 and D = 6, we have 



-1 + V3, (3.33) 



%min,6 ~ (3.34) 

and 

h{Xh,rrvin) = 12>/3, (3.35) 

2 8 

/e^e.mm) = ^32- (3.36) 
Hence, the boundary problem can be solved in D = 5 and D = 6 only if the conditions 

— > ^vrv/3, (3.37) 

- > ^ (3.38) 
Pe 27 

are satisfied. We thus have the same critical effect as in the four-dimensional dS 
space-time. The trapped surface can be formed only if the energy of colliding particles 
is not very high. Assuming that Xq <C 1, we obtain 

Pa « 2 (^) (3.39) 

from 

3.3 Area of the marginally trapped surface below the critical point 

Knowing p , we can calculate the area of the trapped surface: 

Ara P = 2 y v/det^dV, (3.40) 

P<P0 

where g a p is the induced metric on the trapped surface and dV is an elementary 
volume in the (D— 2)-dimensional flat space. Because the corresponding induced 
metric in the four- and five- dimensional space-times has the form (the form is similar 
in higher dimensions) 



/ 1 \ 

1 / 1 \ 1 



9ap -Jr i o ' 9aP ~Jr 



1 
\ 1 / 



(3.41) 
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the explicit expression for the area of the trapped surface is 

A rap = 2 • Vol S D -* I K J_^ dp = 2 • Vol S D -* I —P—dp. (3.42) 



J\f( D ~ 2 )/ 2 J (l + p 2 /2a 2 

o o 



In particular, 



A* = Sn,^^, (3.43) 



2a 2 + p' 

(3,4, 



We have 



^ « Stt-^ « Sttp 2 f 1 - ^ ) ■ (3.10) 

gt + V or 



-4| E - « (3 ' 46) 

in the low-energy limit and 

A 4 Cl = (4 - 2^)^, (3.47) 

= (l - Sv^arctan ^- 1 ^ a 3 (3.48) 

at the critical point. 

4. Concluding Remarks 

We have studied the formation of marginally trapped surfaces in head-on collisions 
of two shock waves in the dS space-time for D < 3. For D < 4, we found the 
critical value of of the shock wave energy dependent on the dS radius, above which 
the trapped surface equation has no solution. This critical behavior is similar to 



that found in [2TJ and is also reminiscent of the behavior encountered in numerical 
simulations of gravitational collapse For D = 3, the critical energy above 

which there is no trapped surface is independent of the cosmological constant. 
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5. Appendices 



A. Geometric view of shock waves in plane coordinates 

A single shock wave in the D-dimensional dS space is shown in Fig. 0.A. The dS space 
is represented as a hyperboloid embedded into the (D+l)-dimensional Minkowski 
space-time. The presented shock wave is located on the intersection of the hyper- 
boloid and the plane x° — x 1 = 0, 

/•y»0 I rf 1 / > ■ / >■ 1 



The coordinates x 2 and x 3 are suppressed in this figure. 

Two shock waves colliding at u = v = are shown in Fig. |2].B. 




Figure 2: A. A single shock wave in the dS space can be represented as the intersection 
of the hyperboloid and the plane x° — x 1 = (x 2 and x 3 are suppressed). B. Two shock 
waves in the dS space. A collision of two shock waves occurs at x° = and corresponds to 
the collision of red and yellow balls. 

We can also make an animation and draw the position of the single shock wave 
at discrete instants. In Fig. [| we see that the shock wave is a nonexpanding one. 




Figure 3: A shock wave in the dS space at different instants of u x time." 
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B. Solution of the geodesies equation 
B.l The (j-model (n-field) approach 

In this section, we derive and solve the null-geodesic equations in the dS space- 
time with a nonexpanding shock-wave. We note that applying the embedding the- 
orems |^0| requires a delicate analysis because of the nonsmoothness of the metric 



with the shock wave. 

We use an analytical approach similar to the a- mo del (n- field) approach and 
start from the Lagrangian 



C 



dr 



dxM{T) -G M N(x(r))^l- A (x M (r)g MN x N (r) - a 2 ) 



dr 



(B.l 



where M, N — 0, . . . , D, g MN is the metric defined the hyperboloid, and Gmn{x{t 
is the metric deformed by the shock wave, 

G M n = 9mn + F8{u)du 2 . 



The Euler equations for this Lagrangian are 



d*X»{T) N 
(JMN — h MN*- KL 



dr 2 



dx K (r) dx L (r) N 

— j j VXg M NX (r) 

dr dr 

x M (r)g MN x N (r)-a 2 



0. 
0. 



Using u5(u) = 0, we obtain 

r MN n 

and equations (B.3| ) and (|B.4|) reduce to 



K - x M , 



x N + Y N KT x K x L + \x 



N 



KL 



0. 



A = — -(x, x + Y KL x K x L ) G . 



(B.2 

(B.3 
(B.4 

(B.5 
(B.6 



Here and hereafter, (x } i/)g = Gmn% y ■ The nonvanishing components of the 



connection T^ K are 



Taking (x, x) g = a 2 and G fiu x tl x l/ = into account, we obtain 



A 



2a 2 



-F + x^^S^ii 2 



(B.7 



(B.8 



Substituting this expression in ( B.6| ) and taking (|B.7|) into account, we obtain 

1 

1 



a 



v - \^5'(u)u 2 - 
F 

i 1 - \f^{u)u 2 



2a 2 
1 

2^' 



-F + x L F^)5{u)u 2 u, (B.9 
-F + x i F t i)5(u)u 2 v, (B.10 
-F + x i F >i )5(u)u 2 x i . (B.ll 
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Noting that the right-hand side of ( |B.9|) vanishes, we obtain u = 0. Taking t — u, 
we obtain 

v - l -F5\u) - F^u)^ = -7^(-F + x l F^)v5(u) } (B.12) 

i l - l -F^(u) = ~^(-F + x j Fj)x l 5(u). (B.13) 

We now use the ansatz (by analogy with the case of the shock wave in the Minkowski 
space-time |33|, [3 



V 


= v + v x u + Q(x 3 )6(u) + R(x J )9(u)u, 


(B.14) 


x i 


= x i0 + x a u + x if 6(u) + Si(x J )9(U)U, 


(B.15) 


i> 


= V! + Q{xi)5{u) + R(x J )6(u), 


(B.16) 


X* 


= xu + Di(xl)8{u) + Si^xDOiu), 


(B.17) 


V 


= Q(4)5'(u)+R(4)5(u), 


(B.18) 


x l 




(B.19) 



Along the geodesies, we have the identity 

= F(xq)5'(u) - (B.20) 

(a similar identity for the flat space-time case was used in fl3~4| ). Transforming our 
anzatz in order to cancel S 2 (u), we now obtain 

v = v + v x u + Q{x j Q )9(u) + R(xi)6(u)u, 
. . . (B.21) 
x % = Xq + x\u + Si0{u)u 



with the bounds 



xl = a 2 , (B.22) 
v = xlx\, (B.23) 

vi = \xf . (B.24) 



And hence 



Q = If, (B.25) 
R = \fA + ^(F- 4*»o + ~Fj + ^(F 2 - KF,) 2 ), (B.26) 
S l = l -F, + ^- 2 {F-x j F,)xl (B.27) 
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We can take x\ = for simplicity. This gives Vq — V\ — 0, and we have 




R = lF 2 t + -±- 2 (F*-(xi ) F t f), (B.28) 

o oCL 

S i = -Fi + -^-(F - xlFAxL 
B.2 Focusing of geodesies 

Before deriving the explicit expression for the trapped surface, we study the structure 
of geodesic beams in terms of the dependent plane coordinates. We consider points 
that are initially on the surface x\\ + x| + x 2 m = a 2 and change to the angular 
coordinates: 



x 20 = a ■ sin ■ sin 6, x 30 = a ■ sin <p • cos 6, x 40 = a ■ cos 0, 
F = AV2p (-2 + ^ln( ? + *!>)) = A V2 P (-2 + cos^ln ( ] + C ° S j )) . 

It is obvious that the term Si(xo)9(u)u in formula ( |B.17| ) leads to the refraction of 
the geodesies. The refraction coefficients are S^Xq). For the refraction coefficients 
in the angular parameterization, we have 

S 2 = -V2^, S 3 = -V2^, S 4 = ^l„fi±i^y (B.29) 
a sin a sin <p 2 a \ 1 — cos / 

The coordinates X2 and £3 correspond to the physical coordinates related to the size 
of the beam. Hence, we can easily find the value of the affine parameter u, which 
determines the focal point of the beam (the point X2 = £3 = 0): 

0,2 • 2 , 
u = ■= ■ sin (p. 

V2p 

The "focal length" is different for each value of the parameter 0, i.e., for each circular 
ring of the beam cross section. This statement is illustrated in Fig. |j. Of course, this 
only roughly explains why the trapped surface can exist. To find a physical meaning 
for this consideration, we must change coordinates to independent ones. 
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Figure 4: Focusing of geodesies: in the dependent coordinates, the focal length changes 
along with the initial conditions. 

C. The shock wave in the independent coordinates 

C.l From plane coordinates to independent coordinates 

To study the structure of the space-time in terms of the independent four- dimensional 
coordinates, it is convenient to use the complex conformal flat coordinates 

2au 2av . \[2a , , . . _ „ . 

™ = ^r—> a = ^—i ( = ^—(x 2 + ix 3 ). (c.i) 

In these coordinates, the shock wave metric is 

2 _ -2dw da + 2d( d( + 2H(£, ()5(w) dw 2 
~ [1 - (w - CC)/2a 2 ] 2 ' 

where 

^■o-K i+ ^ f )'(^S) 



and F is given by (|2.4j). Hence 



,«c,c>^( 1 + ^)(- 2 + i^ ln (f)). «U, 

C.2 An analogue of the D'Eath and Payne coordinates 

(from the coordinates u>, a, ( to the coordinates W, S, T) 

To eliminate 8{u) from the metric, we use a coordinate change analogous to the 



change introduced in f45 



w = W, 

a = Z + H(r i r)6(W) + W6(W)H r H f} (C.5) 
C = T + W9(W)H?, 
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where Hy = dxH(T, T). In these coordinates, we obtain the metric 

, 2 -2dW + 2\dT + W6(W)(H T y<n + H^dT)\ 2 . 
ds = Fi 777777 , 777777777777^ /0 9 19 > ( C - 6 ) 



[1 - (WE - XT + W0(W)G)/2a 2 ] 2 

where G = H — TH T — Tif^ an d H(T, T) depends on T, T as (C, C) given by ( |C.3| ) 
depends on (,(. 

C.3 Geodesies in terms of the independent complex conformal flat coor- 
dinates 

We have the expressions for geodesies in terms of the dependent coordinates. Us- 
ing coordinate change flU.lp , we obtain an expression for geodesies in terms of the 
independent coordinates. In the first order of the parameter u, we obtain 

w(u) = W\U + . . . , 

a(u) = cr + aiu H = a 0c + a g8(u) + {a Xc + a ie 6{u))u + (C.7) 

C(u) = Co + &« + • • ■ = Co c + CoeO(u) + (Cic + CwO(u))u + . . . , 

where 

wi = rrV- ( c - 8 ) 



1 + Xp/a' 
1 + X(j/a' 



^oc = - , ° , (C.9) 



2Q(4 

1 + x%ja 



^oe = — — rr, (C.10J 



<7lc_2 U + ^/a a(l + xVa) 2 ) ' ( °' n) 

o ( g(4M QsVo) £ 4 (^o>o \ t (ai2) 

l+Xg/a a(l + Xo/a) 2 a(l + Xg/a) 2 a(l + XqCl) 2 J 



V2z 
1 + Xq/o' 

Cos = 0, (C.14) 

\f2zi x\ V%z 
1 + x^/a a (I + Xq/o)' 



Coc = t~~ 4 ° , (C13) 



Clc - i I 1471 TTl I ^4/„\2> (C.15) 



Ue 1+^/a a(l + x^/a) 2 ' 1 J 
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where the complex variables S, Zq, z\ are related to S\ Xq, x\ by 

S = S 2 + iS 3 , (C.17) 

zo = xl + ixl, z\ = x\ + ix\. (C.18) 

We see that there is a discontinuity only for the a variable. 

C.4 Geodesies in terms of the independent smooth coordinates 

Using relations (|C.5|) between the initial independent coordinates and the smooth in- 



dependent coordinates, we obtain the expression for geodesies in terms of the smooth 
independent coordinates up to the second order: 

7/1 2 

E(w) = (E 0c + E oe 9(w)) + (E lc + E ie 9(w))w + E 2 — + . . . , (C.19) 

w 2 

T(w) = T + (T^ + T w 6(w))w + T 2 — + . . . , (C.20) 

— — — — — tl) 

t(w) = To + (T Xc + t w 0(w))w + T 2 — + . . . , (C.21) 

where 

S 0c = a , (C.22) 

Z oe = 0, (C.23) 

Si c = (Tic, (C.24) 
E 10 = a w - H r (T , T )T 1 - tf T (T , t )T 1 - w 1 H r (T , f )Hr(T , T ), (C.25) 

T = Co c , (C26) 

Y lc = Cic, (C.27) 

Y w = Ci e -^i^T(T ,T ). (C.28) 

D. Trapped surface equation 

In the case of two shock waves, we deal with 

ds 2 = 2gwzdWdE + gyydTdT + 2g T fdTdT + gyydTdT = 

-2dWdE + 2|dT + (E9(E) + W9(W))(H rr dT + H ff dt)\ 2 



[1 - (WE - TT + (E9(E) + We(W))G/2a?)]< 



(D.l) 



Because of the notation symmetry (W <-> E), we may analyze only one part of the 
trapped surface. 
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The null geodesies passing through the surface 

£ = tf(T,T), W = (D.2) 
can be specified by the tangent vectors (see Subsections C.3 and C.4): 

£ w = Wl , e = °i, e T = Ci, e f = Ci, (d.3) 

where ^ 

o-i = ~o (0ttCi + 2^ T xCiCi + 9fr(i)- (D.4) 

We also suppose that 

(£,#a)=0, a = C,C, (D.5) 

for 

= (0, -0 a ¥, ({) (D.6) 
(our calculations are very close to those in the review section in J28[). This gives 

-0WE0.tfWi + 0a6C 6 = O. (D-7) 

From ( p.7| ), we have 

C = w^gw^, (D.8) 

where 

g ab ig 

the inverse metric for g a b. Therefore, we have 

i W = w u (D.9) 

£ s = o x = Wl3 Y\-, (gyydy^dy^ - 2g rf d T ^df^ + g rr df^df^), (D.10) 

2aet(g ab ) 

£ T = G = j^fpFj (9tfdr^ - grrdyy), (D.ll) 

f = Ci = ^^(^tt^* - ^tt At*)- (D.12) 
For £m = 9mn£, N at the point W = 0, we have 

/ dy^dy^ 1 Oy^ dy^ \ 

V 1 + TT/2a 2 '~l + TT/2a 2 '~l + TT/2a 2 '~l + TT/2aV ' ^ D ' 13 ^ 
For the convergence at the surface W = 0, E = — £), we obtain 

= /i M7V V^a/, (D.14) 

where 



.MAT 



K™g af3 K$ (D.15) 



- 19 - 



and g a P is the inverse of the metric induced on the trapped surface, 



g af 3 = K^QMNKp 



g rr g rr \ 
Qtt QrtJ 



(D.16) 



The components of the connection for metric (p.l|) in the coordinates X M 
(W, E, T, X), M, N = 0, 1, 2, 3, are 



ii 



r 1 

1 o 



22 



r 3 

1 13 



r 1 
1 1 



2J\f ' 

II YT 



r} 



dyAf 



r 2 

1 22 



12 2jV 



r 2 

1 12 



r 1 

1 23 



33 



r 3 

1 33 



AA 



2JV 
<9 T Af 



2A/" ' 



13 



2A/" 



23 



<9 s Af 
' 2AT ' 



AT 



where 



A/"= [l - (WE -TT + (E0(E) + ^(^))G)/2a 2 ] 
At the point W — 0, the tensor /i M7V can be represented as 

/o \ 



h 



MN 



(D.17) 
(D.18) 
(D.19) 
(D.20) 

(D.21) 



(D.22) 



-2d T ^df^ <9x^ <9 T ^ 
det{g al3 ) <9 T ^ 0-1 

\0 <9 T ^ -1 / 

Finally, the equation for the trapped surface can be found by direct calculations: 

2 



1 + 



2^2 



<9"f T "I 9 



2^-# 
1 + TT /2a? 



0. 



(D.23) 
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